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ABSTRACT 


The existence of highly correlated quasi-biennial variations 
in the geomagnetic field and in solar activity Is demonstrated. The 
analysis uses a numerical filter technique applied to monthly averages 
of the geomagnetic horizontal component and of the Zurich relative 
sunspot number. Striking correlations are found between the quasi- 
biennial geomagnetic variations determined from several magnetic 
observatories located at widely different longitudes, indicating a 
worldwide nature of the obtained variation. The correlation coef- 
ficient between the filtered Dst index and the filtered relative 
sunspot number is found to be -0.79 at confidence level >99% with a 
time-lag of 4 months, with solar activity preceding the Dst varia- 
tion. The correlation between the unfiltered data of Dst and of the 
sunspot number is also high with a similar time-lag. Such a time- 
lag has not been discussed in the literature, and a further study 
is required to establish the mode of sun-earth relationship that 
gives this time delay. 
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INTRODUCTION 


Oscillations in various geophysical quantities having periods of 
about two years, which are referred to as quasi-biennial variations 
in this paper, have received considerable attention in recent years. 

In the case of such a variation in the geomagnetic field its exis- 
tence and its exact nature are not well established. This is mainly 
because of the difficulty in extracting long period, small amplitude 
variations from records that contain other changes of larger ampli- 
tudes such as the daily or seasonal variations and disturbances of 
all sizes. The limited availability of magnetic observatory re- 
cords having a long time span and high accuracy is another factor 
for the lack of definitive studies. Nevertheless, the study of 
quasi -periodic variations in the geomagnetic field is of interest 
because it provides valuable information concerning the solar- 
terrestrial relationships on one hand, and about the earth's 
interior on the other. Quasi -periodic variations in meteorological 
parameters have recently been extensively investigated because of 
general interest in long-term climatic changes and in solar weather 
relationships (e.g., King, 1975; Mason, 1976; Landsberg and Kaylor, 
1976). 

The interactions between the sun and the earth's magnetos- 
phere and atmosphere have been studied intensively since in situ 
measurements in space became possible by means of satellites. Be- 
cause of the relatively short time span of this new space exploration 
era these studies have been almost exclusively on short term solar- 
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terrestrial interactions. However, the long term interactions between 
the sun and the earth's environment may be of nature different from 
interactions for short lived disturbances such as magnetospheric 
substorms or magnetic storms. Therefore, it is of interest to in- 
vestigate long term variations, for instance, in the geomagnetic 
field and to seek the causal relationship with solar phenomena. Such 
a study may lead to an understanding of different modes of the solar 
influences on the earth's environment. 

BACKGROUND 

The first mention of a quasi -periodic geomagnetic variation with a 
period of about two years appears to have been made by Kalinin (1952, 
1954). However, Kalinin's results were regarded as being inconclusive 
by Currie (1966) because Currie's power spectral analysis by means of 
the Blackman and Tukey (1958) autocorrelation method showed no spectral 
peak near two years of period. On the other hand, Mansurov (1960) de- 
rived the secular change for the period of 1957 to 1959 for the Soviet 
Antarctic station Mirny by the same analysis procedure as that employ- 
ed by Kalinin, and showed the presence of about one and one-half cycles 
of a quasi-biennial variation. In a more recent power spectral analy- 
sis of geomagnetic data by meins of the maximum entropy method, Currie 
(1973) identified the spectral bands of the double solar cycle (average 
period 21.4 years) and the solar cycle (10.6 years) in the power 
spectral densities in the geomagnetic H and Z components. He inter- 
preted the spectral peak near 2.15 years as the ninth and fifth 
harmonics of the double solar cycle and the solar cycle, respectively. 
Fraser-Smith (1972) gave a negative result on the existence of a quasi - 
biennial variation in the geomagnetic activity index Kp. 
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With a stlmulous from the discovery by Reed and Rogers (1962) of 
the "26-month M zonal wind oscillation In the equatorial stratosphere, 
Shapiro and Ward (1962) made a spectral analysis of the ZUrlch rela- 
tive sunspot number and pointed out the possibility of the existence 
of a spectral peak at the period of 25 months. However, their result 
Is not conclusive since this peak was shown to be significant only at 
5 percent level when subjected to a standard statistical test. In 
their search for similar oscillations in the geomagnetic field, Stacey 
and Westcott (1962) analyzed data on the horizontal component, H, of 
the magnetic field at Huancayo (12° S geographic latitude), Allbag 
(19° N), Apia (14° S), and Eskdalemuir (55° N) by the same method 
as that used by Kalinin (1554), namely by taking the difference be- 
tween the 12-month and 24-month sliding averages of the monthly 
means, and obtained similar results to those derived by Kalinin. 
However, Stacey and Westcott stated that the correlation between the 
three equatorial stations was not consistently good. They performed 
a power spectral analysis, which indicated power spectral peaks near 
26 to 27 months for the Alibag and Apia data; but the corresponding 
power spectral peak for Huancayo was weak and Ill-defined. Pur- 
suing the same question, London and Matsushita (1963) found no 
spectral peak near 26 months in the monthly means of quiet day ampli- 
tude of H from Huancayo. 

Recently, Sugiura (1976) published a preliminary study (hereafter 
referred to as Paper 1), in which he showed evidence for the existence 
of quasi-biennial variations in the geomagnetic field and in solar 
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activity as represented by the Zurich relative sunspot number, and fur- 
ther Indicated a statistically significant correlation between these 
variations. He suggested the possibility that quasi-biennial varia- 
tions observed in the geomagnetic field, cosmic rays, and stratospheric 
and meteorological parameters are produced' by a common cause on the sun. 

In Suglura's analysis, annual averages were used both in the solar and 
geomagnetic data; therefore, the phase and amplitude of the quasi - 
biennial variations were not accurately determined. In the present 
paper the analysis is greatly Improved by using numerically filtered 
monthly averages of solar and geomagnetic data. It is shown that the 
correlation between the quasi-biennial variation in the geomagnetic 
field and that in the relative sunspot number is indeed extremely high 
with a time lag that suggests a different mode of solar- terrestrial 
interaction than those presently known. 

METHOD OF ANALYSIS 

To obtain the amplitudes and phases of the quasi -biennial variations, 
time series consisting of the monthly averages of the relevant solar and 
geomagnetic data were passed through a band-pass filter. This numerical 
filter is the same as that used by Landsberg et al . (1963) for their 
analysis of the surface temperature. Characteristics of the filter are 
briefly described in the Appendix. A comparison of the present tech- 
nique with the method of analysis used in Paper 1 is also made in the 
Appendix. This comparison is of value because both techniques offer 
respective advantages. For instance, while the former method provides 
detailed information on the wave form, the latter technique utilizing 
the second order time derivative of annual averages reduces the required 
numerical analysis by a large factor, and still preserves the coherency 
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between data from different observatories. 

DATA 

The magnetic observatories whose data are analyzed in this paper, 
their coordinates, and the periods for which monthly averages were used 
in the analysis are given In Table 1. The records from Cheltenham and 
Fredericksburg were regarded as a single time series by subtracting a 
constant value from the Fredericksburg data so as to make the long-term 
variations as continuous as possible. Likewise, whenever an observatory 
moved Its location, a constant value Is added or subtracted to form a 
single time series. 

QUASI-BIENNIAL VARIATION IN GEOMAGNETIC FIELD 

Figure 1 presents filtered data of the horizontal component, H, 
of the magnetic field obtained at San Juan, Cheltenham-Fredericksburg, 
Tucson, Honolulu, and Kakloka, clearly demonstrating the presence of 
persistent, coherent quasi-biennial oscillations. The periods of these 
oscillations are about 2 years, but there are considerable variations 
in the time interval between neighboring peaks, ranging from about 1% 
years to 3 years. In some cases, the apparent change in the peak-to- 
peak Interval or an Irregular wave form Is due to a superposition of 
two trains of waves with a phase lag. However, the similarities among 
the oscillations observed at these observatories are, on the whole, 
striking. Indeed, the similarities are more impressive than are shown 
In Paper 1 on the basis of the second time derivatives of annual averages. 

To show the high correlations between the data from pairs of obser- 
vatories, the correlation coefficients were computed, the results of 
which are given in Table 2. In the table, "Student's t and the number 
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of degrees of freedom, n, are provided besides the correlation coefficient 
to enable one to judge statistical significance. On the basis of the 
standard t-test, the probability that the obtained correlation should 
arise by random sampling from an uncorrelated population Is much less than 
0,001 In all cases; see, e.g., Fisher and Yates (1957) for the "Student's" 
t distribution. There can be no question about the presence of coherent 
quasi-biennial oscillations at these observatories. 

In Paper 1, the amplitudes of the quasi-biennial variations derived 
from Kakloka data for disturbed, all, and quiet days were compared 
(Figure 2 In that paper). The amplitudes so obtained were then found to 
diminish in the order: disturbed days, all days, and quiet days. If the 
source of the quasi-biennial variation lies in the interior of the earth, 
one would expect the amplitude to be independent of magnetic activity which 
Is of origin external to the earth. Therefore, It was concluded that the 
observed oscillations are caused by an external source. It was then shown 
that the second time derivative of the Dst index and that of the Zurich 
relative sunspot number were correlated at a statistically significant 
level (Figure 3 in Paper 1). It is of interest to examine this correla- 
tion with numerically filtered monthly averages of these activity Indices 
because of the improved time resolution achieved by this technique. Be- 
fore we make such a comparison we first demonstrate the existence of a 
quasi-biennial variation In solar activity. 

QUASI-BIENNIAL VARIATION IN SOLAR ACTIVITY 

By an application of the band-pass filter to the monthly averages 
of the Zurich sunspot number for the years 1900 to 1976, the existence 
of persistent quasi-biennial oscillations in this solar activity index 
was indicated as clearly as in the geomagnetic field. (Examples are 
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shown in Figures 2 and 3). As In the geomagnetic data the period of 
oscillation and the amplitude vary somewhat, but the regularity of the 
oscillations was striking. When the same filter Is applied to the 
Ottawa 10.7 cm solar flux, we again obtain evidence of a quasi-biennial 
variation. In fact, as can be seen In Figure 2, the quasi-biennial 
oscillations in 10.7 cm solar flux are sc coherent with those In the 
Zurich relative sunspot number that the two curves are almost Identical. 
This result is not suprising because the relative sunspot number and 
10.7 cm solar flux themselves are highly correlated. With the data for 
1947 to 1975 the correlation coefficient between the monthly averages 
of the Zurich relative sunspot number and of the Ottawa 10.7 cm solar 
flux was found to be 0.983 with n=345 and £=100.7, where n and t are the 
number of degrees of freedom and "Student's" t. The correlation coef- 
ficient between the filtered data of these parameters is 0.949 with 
n s 291 and t_=51 . 1 . (The number of pairs is less for this set because of 
the loss of the data points at both ends arising from the application 
of the numerical filter). These correlations are significant at >99.9% 
on the basis of t- tests. 

CORRELATION BETWEEN OSCILLATIONS IN Dst AND SUNSPOT NUMBER 

The magnetic activity index Dst represents the axially symmetric 
component of the geomagnetic field variations primarily caused by the 
sun. It mainly represents the field from the magnetospheric equator- 
ial current system which exists even at magnetically quiet times (Suglura, 
1973; Suglura and Poros, 1973). There are smaller contributions from 
the magnetopause current and the tail current. 
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The correlation between the numerically filtered monthly averages 
of Dst and those of the relative sunspot number Is found to be even 
better than the correlation between annual averages of those parameters 
shown In Paper 1. Plots of the filtered Dst and sunspot data are pre- 
sented In Figure 3, which shows a remarkably good anticorrelatlon; 
namely, the phases of high (or low) solar activity correspond to those 
of low (or high) Dst level. It was further found by visual inspection 
that this correlation Is better when the Dst curve Is advanced In time 
by about 4 months, l.e., when Dst Is correlated with the sunspot number 
4 months earlier. To study this phase relation quantitatively we pre- 
sent In Figure 4 the correlation coefficient between Dst and the 
relative sunspot number, R, and "Student’s” t as functions of lag time. 

The maximum absolute value of the correlation coefficient occurs at a 
time-lag of about 4 months, the correlation coefficient being -0.79. 

The value of t is a maximum at this time-lag, meaning that the proba- 
bility that the obtained correlation might occur by a random sampling 
from an uncorrelated population is a minimum for this time-lag. For 
the 4-month time-lag, such a probability is less than 0.001. (The de- 
grees of freedom for the plotted points are between 158 and 172). Hence 
the result obtained here can be regarded as being statistically signifi- 
cant. However, the value of 4 months was determined statistically 
involving about eight cycles of oscillation and the phase differences 
for individual cycles fluctuate slightly, but the plots in Figure 3 and 
the results shown in Figure 4 indicate that the time-lag determined here 
is not a mere average of widely varying values, but represents a meaningful 
average value. 

In Figure 4, we further observe that the next maximum in the correla- 
tion comes near a lag of +8 months; the correlation coefficient is +0.75, 
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again coinciding with a peak In t. The time interval of 12 months 
between the two maximums represents one half of the period, confirming 
that the oscillations we are dealing with has a period of about 2 years. 

In interpreting the 4-month time delay It is important to know 
whether this time-lag is peculiar to the quasi-biennial variation or 
the same time-lag exists in the correlation between the original, un- 
filtered sunspot data and Dst. Thus the correlation coefficient was 
calculated between the unfiltered monthly averages of the Zurich rela- 
tive sunspot number and those of Dst with varying time-lag, ?, at intervals 
of 1 month from T= -12 months to T= +7 months and then at b-month 
intervals outside this range. The results are shown in Figure 5, where 
$ now in units of years instead of months as in Figure 4. The maxi- 
mum correlation is found at T* -3 months, where t is also a maximum. 

The number of degrees of freedom is 22b in all cases, and for the region 
of T giving high correlations the confidence level is well above 99%. The 
high (positive) correlation near 'C* -b years is obviously from the out- 
of-phase correlation of the (negative) correlation for the 11-year solar 
cycle change. The correlation coefficient varies slowly from T= -3 to 
-4 months; hence it is not certain whether or not the difference betwe^ 
the 3-month time-lag in this correlation and the 4-month time-lag in 
the correlation of the filtered data is significant. However, it is 
likely thaw the cause for the 4-month time-lag in the quasi-biennial 
variation in Dst relative to that in the sunspot number is the same as 
the cause by which the variations in the monthly average Dst 1-g behind 
those in the monthly average relative sunspot number. 

The results in Figure 5 are somewhat unexpected because it is known 
that the annual mean intensity of magnetic storms peaks 2 or 3 years 
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after the sunspot maximum (e.g., Sugiura and Chapman, 1960). The 
present results mean that statistically, the Integrated ‘Dst effect* 
lags behind solar activity represented by the relative sunspot number 
by 3 or 4 months. A similar result has also been obtained by using 
the geomagnetic activity index Ap instead of Dst. This subject, which 
is outside the scope of this study, is now being investigated in more 
detail . 

DISCUSSIONS 

1. In Figure 1 we have plotted the San Juan data only after 1930, 
though the filtered data for this observatory begins from April 1927. 

The omission of data for the first 2 3/4 years was made because of 
the presence of waves of unusually large amplitude reaching values 
several times the normal amplitudes. We believe that these large 
amplitude oscillations are from artificial sources that existed in the 
early phase of the operation after the observatory moved from Vieques 
Island to San Juan; see the footnote in Table 1. Filtered data were 
obtained for Vieques Island, but were not included in Figure 1 because 
there was nothing unusual in the short stretch of data from this sta- 
tion. While preparing Paper 1 it was noticed that the second order 
time derivatives of annual mean H values indicated extremely large 
amplitude oscillations for some observatories (not shown in Paper 1), 
especially for the early years, say before 1940. It appears that the 
deduction of quasi-biennial oscillations provides a convenient means 
of testing the quality of observatory magnetic data. The application 
of the analysis method of Paper 1 to the vertical component, Z, indi- 
cated that the measurements by the conventional variometers are generally 
much less accurate for Z than for H and that the accuracies in the Z 
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data from many observatories seemed inadequate for the study of the 
quasi -biennial variations. 

2. The oscillations seen in the numerically filtered data must be 
viewed with caution because of the application of the narrow band-pass 
filter (described in the Appendix). However, the existence of geomag- 
netic field oscillations having periods of a few years appears to be 
real. This is because the analyses made by Kalinin (1954) and Sugiura 
(1976) and in this paper utilizing entirely different techniques have 
all indicated the presence of such oscillations. In addition, the 
existence of a spectral peak or band near 2 years in the power spectral 
density distributions derived from various sets of geomagnetic data 
supports the above view (e.g., Currie, 1973). Our own power spectral 
. ilysis of the geomagnetic and solar data used in this paper by means 
of the maximum entropy method (e.g., Ulrych and Bishop, 1975) and also 
by the method of Blackman and Tukey (1958) indicates the existence of 
a band in the period range roughly 2 to 3 years. Figure 6 shows an 
example of such results; namely in this case, the power spectral den- 
sity in the relative sunspot number calculated both by the maximum 
entropy method and by the Blackman and Tukey method. The frequency 
range for which the filter transfer function is greater than one half 
the maximum is indicated by a small box drawn in the figure; crudely 
speaking, oscillations having frequencies in this range are passed by 
the filter used in this analysis. Figure 6 clearly shows the presence 
of power in a band near 2 years. A detailed study of the power spec- 
tral densities in solar and geomagnetic activities is outside the scope 
of this paper and will be discussed in a separate report. 

The correlation with high statistical significance between the 

relative sunspot number and the geomagnetic H component (or Dst) cannot 
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be questioned. However, the exact nature of the signals transmitted 
from the sun to the earth and giving a high correlation is not as yet 
clear. A visual comparison of the unfiltered sunspot and geomagnetic 
data does not offer any clue on this point, and at the present stage 
of the study the filtered data shown in this paper seem to provide 
a good means of presenting highly correlated solar and geomagnetic 
signals. The higher maximum correlation of (-)0.79 between the 
filtered solar and Dst data compared with the maximum correlation 
coefficient of (-)0.52 for the unfiltered data indicates a significance 
of the quasi-biennial variations. 
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CONCLUSIONS 


By means of a numerical filter technique we have shown the exist- 
ence of highly correlated quasi-biennial variations in solar activity 
and In the geomagnetic field. The following conclusions have been 
drawn from the present analysis: 

(i) The quasi -biennial variations in the H component determined 
from several observatories at different longitudes, ranging from 140°E 
to 294°E, are found to have a strong linear dependence among themselves. 

(il) The quasi -biennial variations deduced from the Zurich rela- 
tive sunspot number and the Ottawa 10.7 cm solar flux are virtually 
identical . 

(iii) T‘'e**e is a strong linear dependence between the quasi- 
biennial variati n in the relative sunspot number (or 10.7 cm flux) 
and that in the geomagnetic Dst index which represents variations in 
the axially symmetric external field. The correlation between these 
quantities is a maximum at a time-lag of 4 months with solar activity 
preceding the Dst variation. 

(iv) A similar time-lag is obtained in the correlation between 
the unfiltered data for the Dst and for the sunspot number. 

(v) The origin of the 4-month time-lag in these correlations is 
not known. Since the relative sunspot number is used merely as a para- 
meter expressing solar activity, this time-lag may not represent the 
true time delay in the causal relationship. 

(vi) Estimates of the power spectral density of the relative 
sunspot number indicate the presence of a band of signals near 2 
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years period. However, the exact nature of the highly correlated solar 
and geomagnetic oscillations is not as yet understood. 

Ac know! edgements . We wish to thank Dr. H. E. Landsberg, Dr. J. W. 

King, and Dr. T. Yukutake for valuable discussions and Dr. J. London 
for calling our attention to his and Dr. Matsushita's earlier study. 

We are indebted to W. B. Schar for his assistance in the data analysis. 
We also wish to express our thanks to: Dr. K. Yanagihara and Dr. M. 
Kawamura, respectively the former and present Director of the Kakioka 
Magnetic Observatory for providing the year books and additional data, 
and the World Data Center-A for Solar-Terrestrial Physics in Boulder 
for supplying the geomagnetic and solar data used in this analysis. 


14 



REFERENCES 


Blackman, R. B., and J. W. Tukey, The Measurement of Power Spectra , Dover 
Publications, Inc., New York, NY, 1958. 

Currie, R. G., The geomagnetic spectrum - 40 days to 5.5 years, J. Geo- 
phys . Res . , 71_, 4579, 1966. 

Currie, R. G., Geomagnetic line spectra - 2 to 70 years, Astrophys . 

Space Sci. , 21_, 425, 1973. 

Fisher, R. A., and F. Yates, Statistical Tables for Biological, Agricul- 
tural, and Medical Research , Oliver and Boyd, London, 1957. 

Fraser-Smith, A. C., Spectrum of the geomagnetic activity index Ap, J^ 
Geophys. Res. , 77 _, 4209, 1972. 

Kalinin, Yu. D., Some questions in the study of geomagnetic secular 
variations, in Russian, NIIZM, 8(18), 5, 1952. 

Kalinin, Yu. D., The long-period geomagnetic variation, in Handbook of 
the variable geomagnetic field in the USSR, edited by V. I. Afanas'yeva, 
Leningrad, Hidrometeoizdat, in Russian, pp. 28-43, 1954; Defence Res. 
Board of Canada , Translation T379R , 1954. 

King, J. W., Sun-weather relationships, Astronautics and Aeronautics , 13. 
(4), 10, 1975. 

Lanczos, C., Applied Analysis , p. 305, Prentice-Hall, Inc., Englewood 
Cliff, New Jersey, 1964. 

Landsberg, H. E., and R. E. Kaylor, Spectral analysis of long meteoro- 
logical series, Inst. J. Cycle Res. , 7 _ , 1976 (in press). 

Landsberg, H. E., J. M. Mitchell, Jr., H. L. Crutcher, F. T. Quinlan, 
Surface signs of the biennial atmospheric pulse, Mon. Weather Rev. , 91 , 
Nos. 10-12, 549, 1963. 

London, J., and S. Matsushita, Periodicities of the geomagnetic variation 
field at Huancayo, Peru, Nature , 198 , 374, 1963. 


15 


Mansurov, S. M. , Geomagnetic secular variation in East Antarctica, in 
Geomagnetic Disturbance , in Russian, Moscow, Acad. Sci. Press, pp. 48- 
52, 1960; Defence Research Board of Canada translation T375R , 1960. 

Mason, B. J., Towards the understanding and prediction of climatic varia- 
tions, Quart. J. R. Met. Soc. , 102 , 473, 1976. 

Reed, R. J., and D. G. Rogers, The circulation of the tropical stratosphere 
in the years 1954-1960, J. Atmos. Sci. , 19 , 127, 1962. 

Shapiro, R. , and F. Ward, A neglected cycle in sunspot numbers?, J. Atmos. 
Sci . , 19., 506, 1962. 

Stacey, F. D., and P. Westcott, Possibility of a 26- or 27-month periodi- 
city in the equatorial geomagnetic field and its correlation with strato- 
spheric winds, Nature , 198 , 730, 1962. 

Sugiura, M., Quiet time magnetospheric field depression at 2.3 - 3.6 R £ , 

J. Geophys. Res. , 78 , 3182, 1973. 

Sugiura, M. , Quasi-biennial geomagnetic variation caused by the sun, 
Geophys. Res. Letters , 3^, 643, 1976. 

Sugiura, M. , and S. Chapman, The average morphology of geomagnetic storms 
with sudden commencement, Abhandl. Akad. Wiss. Gottingen Math. Phys. K1 . , 
Sonderheft Nr. 4, Gdttingen, 1960. 

Sugiura, M., and D. J. Poros, A magnetospheric field model Incorporating 
the Ogo 3 and 5 magnetic field observations, Planet. Space Sci. , 8, 921, 
1973. 

Ulrych, T. J., and T. N. Bishop, Maximum entropy spectral analysis and 
auroregressive decomposition, Rev. Geophys. Space Phys. , 13 , 183, 1975. 

Yukutake, T., The solar cycle contribution to the secular change in the 
geomagnetic field, J. Geomagn. Geoelec. , V7, 287, 1965. 


16 





Table 1. Magnetic observatories, their coordinates, and the years for which monthly average data 
were used. 
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Table 2. The correlation coefficient, f , "Student's" t, and degrees 
of freedom, n, for each pair of observatories; the data 
used consist of filtered monthly averages of the H component. 
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APPENDIX: Numerical Filters Used in This and 
Previous Analyses 

The numerical filter used is a band-pass filter with 55 weights. 

As can be seen In Its transfer function shown in Figure 7, the filter 
very effectively passes only those signals having periods near 2 years. 
The transfer function, Y, Is a maximum at 25 months in period, and is 
reduced to one half the peak value at about 37 and 18 months; at 50 
and 15 months in period, Y is down to one- tenth the peak. The side 
lobes are less than 3.4% of the maximum Y. Since the weights are sym- 
metric, no phase shift is Introduced. The numerical values of the 
weights are given in Landsberg et al. (1963). 

In Paper 1 we obtained the second order time derivatives of annual 
averages of the solar and geomagnetic data. Because of the space limita- 
tion in that communication, the details of the analysis technique were 
not discussed. Therefore, we give a brief account of the method and 
compare It with the technique used In this analysis. The previous 
method also can be considered as being a numerical filter technique. 

The underlying Idea is (a) to take annual averages to filter high fre- 

quency components and (b) to amplify signals by a factor w near the 

2 

period of 2 years. The latter comes from the factor co in the second 
time derivative of a wave having an a> ular frequency of <a. With the 
Nyquist frequency at 0.5 cycles per year, or 2 years in period, the tech- 
nique can potentially be useful in extracting quasi-biennial variations. 
However, the effectiveness of the technique depends on the method of 
numerical differentiation. A numerical differentiation is usually made 
by a least squares fit of data to a parabola of second order (e.g., 
Lanczos, 1964). For instance, Yukutake (1965) applied this technique 

to his study of the solar cycle variation in the geomagnetic secular 
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change. However, since the least squares fit, which could be mad -2 with 
any number of points above 3, smooths the data, the method Is not suit- 
able for the present purpose. Therefore, In Paper 1 an exact fit was 
made to a parabola, which was then differentiated twice. This proce- 
dure is equivalent to the simple 3-point difference scheme for the 
second derivative. The transfer function of the procedure is given 
in Figure 8. The figure shows that the filtering of low frequency sig- 
nals Is not very efficient. For comparison, the normalized transfer 
function,(^/ir) 2 * for a heuristic estimate Is also shown In Figure 8 by 
broken lines. The slightly lower effectiveness of the actual numerical 
procedure compared with the ideal case is due to the use of annual 
averages which obviously does not provide adequate time resolution. 

This shortcoming counter-balances the filtering of the unwanted high 
frequency noise by taking annual averages. 

When the results such as those shown in Figure 1 obtained by the 
present technique are compared with the more corse data derived by the 
method of Paper 1 using annual averages, we find that the phase infor- 
mation is lost considerably by the latter method as can be expected 
from the nearness of the Nyquist frequency to the frequencies of the 
signals in question. However, the filter of Paper 1 affected the data 
from different stations and for the sunspot number in the same manner so 
that the correlation between them was preserved to a high degree. None- 
theless, for the study of quasi-biennial variations the use of annual 
means is not recommended except for a crude survey of the presence or 
absence of the oscillations without requiring detailed information on 
the waves. 
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Figure 1. Filtered H date from San Juan. Cheltenham- Frederick sburg. Tucson, Honolulu, and Kakioka, indicating strongly 
correlated quasi-biennial variations. Tide-marks along the vertical axes represent 5 gamma scales far H. 



OTTAWA 10-7 CM SOLAR FLUX 
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Figure 2. Filtered date for the Ottawa 10.7 cm solar flux and for the ZCkich relative sunspot number. 








FiflorB 3. Filtered data for D*t and for tht Zurich relative sunspot number. The Dst 
plots are shifted to the left by 4 months. 
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Figure 4. The coefficient, p, of correlation between the filtered Dst and the filtered relative sunspot 
number and the "Student's" t, as functions of the time-lag, r. 
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Figure 5. The coefficient, p, of correlation between the unfiltered Dst and the unfiltered 
relative sunspot number and the "Student's" t, as functions of the time-lag, t. 




POWER SPECTRAL DENSITY 



Figure 6. The power spectral densities in the ZUrich relative sunspot num- 
ber, estimated by the maximum entropy method and by the Balckman and 
Tukey autocorrelation method. Periods (in years) corresponding to several 
spectral peaks are indicated. The width of the half maximum of the band- 
pass filter is shown by the rectangular window. 
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